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ABSTRACT 


We have investigated the temporal behavior of the x-ray emission 
from Sco X-l on 3 March 1969, based upon data obtained over a continuous 
122 sec exposure in tha energy range 2-20 keV, where the count was 
recorded each 322 jusec. We find no evidence for any temporal structure 
in the data. The 5 o upper limits on the pulsed fraction at any given 
period ?.re below 1% for periods in the range 3 - 300 msec. We also 
establish upper limits for non-periodic variations over times from a 
fraction of a msec to about 1 minute. Any continuous x-ray flickering on 
the scale of seconds is less than 1,27 0 , which is below the 0.02 mag. 
flickering level previously observed in the visible. 


I. INTRODUCTION 


Temporal variations in the intensity of a cosmic x-ray source 
provide critical evidence regarding its physical nature. The outstand- 
ing such example is th^ pulsar in the Crab Nebula, NP 0531, which has 
greatly advanced our understanding of the energy balance for that 
nebular supernova remnant. Another source, Sco X-l, has exhibited 
many variable features in the optical, x-ray and radio (Hiltner and 
Mook, 1967; Lewin, et al. 1968; Abies, 1969). 

The best studied intensity fluctuations from Sco X-l have been 
In the visible. One class of such fluctuations consists of flare- 
type increases with typical amplitudes of about 0.2 mag and with time 
constants of about 10 minutes (Hiltner and Mook, 1967; Westphal et al. * 
1968). Sco X-l was observed to be in a flaring state about 127. of 
the time. In addition, Westphal et al. (1968) found this source to 
flicker, almost continuously, on time scales from several seconds to 
minutes with amplitudes averaging about 0.02 mag. 

Slower variations in the x-ray intensity from Sco X-l are now well 
established. A series of simultaneous measurements from rockets (X-ray) 
and ground-based telescopes (visible) has indicated, on a statistical 
basis, the existence of some correlation in the corresponding gross 
intensities (Chodil et al., 1968; Mark et al. , 1969). During a balloon 
flight, Lewin, et al. (1968) have observed a four- fold increase in hard 
x-rays from this source, with a temporal evolution not inconsistent with 
that of optical flares. More recently, Hudson et al. (1970) have used 
x-ray data from an Orbiting Solar Observatory to correlate two such x-ray 


events with optical counterparts. They found that a two- fold increase 
in the x-ray intensity corresponded to a 0.2 mug optical flare. 

The optical identification of Sco X-l was accomplished via the 
assumption that both the x-ray and optical radiation fluxes arise from 
bremsstrahlung emission by the same isothermal hot, plasma (Sandage, et 
al., 1966), at least under quiescent conditions. Such correlation 
between the x-ray and optical intensities during flares would indicate 
that they are intimately connected even under conditions far from 
equilibrium. Hudson et al. (1970) have suggested that the flare emission 
in both bands is the result of bremsstrahlung radiation ’yy the same 
plasma, where the relative intensities as observed are altered via the 
attenuation of the visible output by self-absorption at the source and 
interstellar extinction. If the x-ray intensity tracks the optical for 
the more rapid fluctuations (flickering), with relative amplitudes 
comparable to what is seen in flares, we can then expect to find x-ray 
flickering amounting to about 10 - 20% of the average intensity. 

However, Westphal et al. (1968) suggest that the optical flickering is 
probably connected with the energy "pump" heating the plasma and that 
it may not show detailed correlations with the x-rays from the plasma. 

On the other hand, the characteristic time scale for energy loss by the 
plasma might be much smaller than that of the optical flickering 
measured to date; Neugebauer f “ al. (1969) find large absorption at 
infrared wavelengths which, if the isothermal bremsstrahlung hypothesis 
is correct, would indicate that the radiating region of plasma is only 
~ lO^cm in size at a density of ~ lO^cm"^, thereby providing a relaxa- 
tion time for energy loss that is < 100 ms. This should certainly 
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give rise to observable x-ray fluctuations if the characteristic 
time for energy input is as large as the time scale for the observed 
optical flickering. We might also expect to see observable short- 
term fluctuations if there is a pulsar in Sco X-l, as has been 
suggested by Tucker (1969). 

We have investigated the temporal behavior of x-rays from Sco X-l 
detected on 3 March 1969 during a continuous 122 sec exposure to this 
source centered at 1003 UT. The results of a spectral analysis cf the 
same data and a description of the experiment have been given else- 
where (Holt, et al., 1969, Boldt, et al., 1969). The temporal informa- 
tion is in the form of the contents of a 3-bit scaler with accumulation 
and readout times of 302 and 20 msec, respectively. The detector count 
rate averaged about 3000 events/sec while pointing at the source. 

During the exposure, the angle between the detector axis and Sco X-l 
was always well within 1 degree, more than an order of magnitude smaller 
than the angular response of the detector. 


II. PROCEDURE 


For a data stream for which the mean value of the instantaneous 
number of counts per temporal bin is N , we can define: 
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If the data stream is temporally structureless, i.e., if the fluctua- 
tions in N* can be attributed solely to Poisson statistics, the variance 

can be formally set equal to N . Furthermore, if the number of 
_ 

temporal bins r is large (£ 40), the p* distribution is approximately 
normal with expectation value (P - 1) and variance 

Let us suppose that the instantaneous counting rate is the 
superposition of a base rate and a disturbance in the form of an 

independent Poisson sample which is present for only some fraction p 
of the time: 


K- -v . ( 2 ) 

The mean M is then: 

H = + O) 

where is the mean which would have been observed had there been 
no disturbance, and 4 may be interpreted as the consequent fractional 
increase in the unperturbed average source intensity. 

Equation (1) may now be evaluated with the expression for M in 
equation (3), and the hypothesis of tested for consistency. 

Should but , the increase in above the expectation value 



An expression for the statistical significance of such a disturbance 
may be stated In t.erms of the number of standard deviations i\. from, 
the ^•■*'0 expectation value for 9(. an< * the total number of counts 
observed ~ • 

/ £ V" /JL_\ Ay/zC 1 

OTTJ 0-f) <V • (5) 

Equation (5) may be interpreted as follows: considering first the 

dependence on |^> , we see that a smaller pulsed fraction can be detected 
by this technique if the duty cycle of the disturbance is small (viz. 
the disturbances are sharp spikes). Similarly, the sensitivity to 
fluctuations increases as the total number of counts (^becomes larger. 
Finally, equation (5) states that, given a total number of events, a 
higher sensitivity to fluctuations can be achieved by lumping the data 
in fewer bins (provided that the bin size does not exceed the on-time 
of the disturbance). Note that if the condition that the bin size 
chosen is smaller than the on-time is satisfied, the sensitivity is 
independent of whether or not the disturbance is contiguous or composed 
of many disturbances such that the net is the same. This is particu- 
larly useful when searching for longer-term fluctuations, and it also 
applies to the important case of periodic fluctuations where coherent 
superposition of data according to a known or a test period results in 
a reduction in the number of temporal bins by the number of folds in the 


data stream. 
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III. NON-PERIODIC VARIATIONS 

Considering first non-periodic fluctuations, we have examined 
the data on the basis of temporal units ranging from 20 msec to 2.6 
sec. A finer temporal division, although possible, would have resulted 
in an insufficient average number of events per bin for the direct 
application of the technique described above. The largest unit of 
time is dictated by the temporal modulation of the aperture provided 
by a paddle in front of the detector; this paddle was advanced to a 
new position every 2.6 seconds. There were five such distinguishable 
positions, each correq> onding to a different effective area. Each 
paddle position was repeated about 10 times during the exposure to 
Sco X-l. In Figure 1, we have plotted the 2.6 sec. rates normalized 
in order to remove the average rate modulation caused by the paddle. 
Standard errors at each point average somewhat more than 1 percent. 

Also shown in Figure 1 are the measured angles between the experiment 
axis and Sco X-l; errors in this measurement have been estimated at 
+ .15°. As the initial gradual increase in the counting rate is caused 
by atmospheric absorption early in the flight, we have included ? n this 
analysis only data to the right of the arrow in this figure. The total 
number of events in the remaining 122 seconds are 333,879; the number 
of 2.6 sec units is 47. We have concluded that any deviations frcm 
constant intensity on this time scale may be attributed to pointing 
errors alone. Using equation (5) with n = 5, we may set a 5 a upper 
limit on the non-periodic pulsed fraction as 
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which, for 0 ■ .5 (the least sensitive case), becomes f < 0.012. 

The next smaller time unit has been chosen to be about 180 msec, 
resulting in 14 times as many temporal bins as were used in the 
previous case. The 5 c upper limit here is f ^ .023 for 0 ■ .5. 

The smallest time scale considered with this procedure is with 5922 
20-msec time units. Here the 5 a upper limit on the pulsed fraction 
is f ;$ .04, where we have again taken g ■ .5. We note that 0 ■ .5 
is truly a worst case, a* the method is as sensitive to intensity 
decreases as it is to increases so that 0 need never be > .5 for 
f « 1. Note also that 0 need not be contiguous. 
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IV. PERIODIC F IXATIONS 

In searching for pulsed emission, we huve chosen to use a fast- fold 
algorithm rather than a harmonic analysis because the former allows for 
a more straight-forward determination of upper limits using the general- 
ized procedure which we have thus far applied to non-perlodlc fluctuations. 
We have searched for possible periodicity In the data by fast-folding at 
all periods In the range 3 - 300 msec, and testing each light curve for 
deviations from randomness. We chose not to attempt the analysis for 
smaller periods, as 10% resolution In the light curve Is as poor as we 
felt was reasonable for making meaningful upper limit statements. For 
trial periods in excess of ~ 10 msec, the temporal resolution used was 
never worse than 2.5%. For the entire range of analysis, 3 - 300 msec, 
the 5 a upper limit for the pulsed fraction f at any given period is 
about .01 for g ■ .5. We note that this value is comparable with the 
upper limits on pulsation set by Friedman et al. (1969) using a harmonic 
analysis, but it covers smaller periods than the 70 msec minimum period 
that they considered and includes the 6 ms period predicted by Tucker 
(1969) for Sco X-l. In addition, we again note that the sensitivity of 
this procedure increases for smaller values of 0, in contrast to the 
harmonic analysis approach. 
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V. RAPID FLICKERING 


The finest time resolution possible in our data is set by the 302 Msec 
accumulation time of the scaler associated with each telemetry word. The 
average number of events per readout is about 1 count (with a scaler 


capacity of 7) so that, when searching for fluctuations on scales from 


fractions of a msec to several msec, we find* It necessary to use a 
different type of analysis. This consists of do termifiing the distribution 
of a large number of scaler readouts and the extent to which this distrib- 
ution differs (if at all) from the Poisson distribution, which may be 
written as ^ 

^ , (7) 
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where ^ is the probability of detecting n counts in t sec when the 

average event rate is \ e *ents/ sec. Since the Poisson distribution '.s 
a one parameter distribution (viz., average event rate), we are essen- 
tially testing the hypothesis that the data in a certain interval of time 
are described by a constant event rate. The actual situation is more 
complicated because the Poisson distribution applies to the ideal case 
of no dead time. A small but non-negligible dead time accompanies every 
event, and this- dead time alters the probabilities for multiple counts 
during a given scaler accumulation period. We have obtained a modified 
distribution which takes into consideration this dead time and may there- 
fore be used in the analysis described above. The suitable modification 
of the Poisson distribution is 
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where p Is the dead time of the apparatus associated with each detected 
event, in this case, o was measured to be 2.8 usee. 

The y2 values obtained indicate agreement In all cases, which would 

appear to rule out large fluctuations on time scales above 300 use c. 

To obtain a quantitative interpretation of this agreement, we refe*. 

to the simplified "on-off" model introduced earlier for candidate dis- 

2 

turbances. For such a small number of degrees of freedom the x 
distribution is quite asymmetric, so that upper limits for this case are 
more properly stated in terms of confidence levels than in terms of a 
number of standard deviations. In figure 2 we have plotted paired values 
of g and f/S that satisfy the condition that the 98% confidence limit for 
the hypothesis of Poisson statistics is exceeded. Any fluctuation 
characterized by parameters which correspond to a point to the right and 
above the computed curve would have been detected on the basis of this 
technique. A detectable fluctuation, for example, involves a pulsed 
fraction f = .015 with @ * .03, independent of how this "on-time" is 
temporally distributed. 
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VI. CONCLUSIONS 

In summary, we have measured the x-ray intensity of Sco X-l and 
found no evidence for fluctuations, periodic >r not. Any continuous 
flickering is less than about .01 magnitudes on the scale of seconds to 
less than .04 magnitudes on the scale of tens of milliseconds. Shorter 
time scale fluctuations are absent to within the limits given in Fig. 2, 
down to a time scale of about 300 jusec. 

These limits severely constrain energy input mechanisms for Sco X-l, 
In particular, we can conclude not only that there are no x-ray intensity 
variations comparable to the optical flickering on a time scale of seconds 
or greater, but that we observe no plasma relaxation on a time scale 
corresponding to the thermal bremsstrahlung parameters of Neugebauer, 
et al. If Sco X-l is a thermal source, therefore, it would appear that 
the energy input to the plasma must be essentially continuous compared 
to the time scale for optical variations. 

With regard to periodic fluctuations, the 5 c upper limits on pulsed 
emission at any given period, in the range 3 - 300 msec, are about .01 
for a roughly sinusoidal pulse profile in time, and considerably less 
(approximately as yfj-p ) for sharper pulse profiles. 
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FIGURE CAPTIONS 

1. Normalized event rates per 302 /usee, averaged over each 2*638 sec* 
paddle dwell time, as function of time. Also plotted Is the angular 
offset of the detector from Sco X- 1 during this exposure, 

2* 98% confidence level limit for rapid (> 300 jusec) fluctuations of 

duty cycle (3 and pulsed fraction f. For time scales £ 10 msec, the 
procedure outlined in Sqction II gives a more sensitive measure 


for fluctuations. 
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